Spintronic Mechanics of a Magnetic Nanoshuttle 
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We investigate theoretically the prospects for using a magnetic nanoelectromechanical single- 
electron tunneling (NEM-SET) device as an electronic spin filter. We find that strong magnetic 
exchange forces on the net spin of the mobile central dot of the NEM-SET structure lead to spin- 
dependent mechanical displacements ( "spin polarons" ) , which give rise to vastly different tunnelling 
probabilities for electrons of different spin. The resulting spin polarization of the current can be 
controlled by bias and gate voltages and be very close to 100 % at voltages and temperatures below 
a characteristic correlation energy set by the sum of the polaronic and Coulomb blockade energies. 

PACS numbers: 85.75.-d, 85.85.+J, 75.76.+j 
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The ability to generate spin-polarized currents is an es- 
sential prerequisite for spintronics applications^. Ac- 
cordingly, many studies of possible spin-filters capable of 
producing a spin-polarized tunnel current between fer- 
romagnetic or normal metal electrodes have been re- 
ported, involving e.g. thin layers of ferromagnetic^ or 
ferrimagnetios insulators, semiconductor quantum dots^, 
Ceo molecules^, and carbon nanotubes^. Quite large spin 
polarization factors have been achieved at low tempera- 
tures, 44% as in Ref. d and 33% as in Ref. [H being typ- 
ical values at 10 K. However, the degree of polarization 
is much lower at room temperature and the pursuit of 
better spin filters continues. 

In this context it is interesting to consider nanoelec- 
tromechanical spin filters. The discrete nature of the 
electronic charge leads to a strong coupling of mechani- 
cal and electronic degrees of freedom in nanoscale single- 
electron tunneling (SET) structures with movable parts, 
such as a SET transistor with a flexible central island. 
The resulting nanoelectromcchanics is rich in mcsoscopic 
phenomena, affecting both the electronic and the me- 
chanical subsystems (see Ref. for a review). 

A conventional approach to the implementation of spin 
controlled nanoelectromechanics is based on spin depen- 
dent tunneling between magnetic conductors. The prob- 
ability for electrons to tunnel is different in the two spin 
channels simply because the electron density of states in a 
magnetic metal is spin dependent. Hence the tunnel cur- 
rent may be spin-polarized with a number of interesting 
consequences for the electromechanics of, e.g., magnetic 
"shuttle" structures^. 

Qualitatively different phenomena may occur in 
nanometer-sized tunnel structures, where short-range 
magnetic exchange forces can be comparable in strength 
to the long-range electrostatic forces between charged 
elements of the device. There is ample evidence that 
the exchange field can be several tesla a few nm from 
the surface of a ferromagnet4~— and the exponential 
decay of the field means that the corresponding force 
on a single electron spin can be very large. These 



spin-dependent exchange forces can give rise to various 
"spintro-mechanical" phenomena and will be of central 
importance here, where we will show that it may result 
in a nanoelectromechanical spin filter with a theoretical 
efficiency close to 100%. 
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FIG. 1: A movable quantum dot (circles) in a magnetic shut- 
tle device can be displaced in response to two types of forces: 
(a) a long-range electrostatic force causing an electromechan- 
ical response if the dot has a net charge and (b) a short-range 
magnetic exchange force leading to a "spintromechanical" re- 
sponse if the dot has a net magnetization (spin). The direc- 
tion of the forces and displacements (arrows) depends on the 
relative signs of the charge and the magnetization, respec- 
tively. 

The theory will be developed for a generic nanoelec- 
tromechanical "shuttle" device in the form of a SET tran- 
sistor with a central island that is movable along the line 
between a source- and a drain electrode as indicated in 
Fig. [T] (where only one electrode is shown). However, the 
theory applies equally well to, e.g., an extended island 
in he form of a suspended carbon nanotube in tunneling 
contact with two fixed electrodes and probed by an STM 
tip and this is the set-up we will use for our quantitative 
analysis. We assume that the island is a "quantum dot" 
in the sense that spatial quantization only allows a single 
electron level to be populated at small bias voltages^. 
In this case both the electronic energy level on the island 
and the probability for electrons to tunnel to and from 
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the bulk electrodes are affected by a mechanical displace- 
ment of the island. An electron that tunnels onto the dot 
changes both the charge and (by its spin) the total mag- 
netization of the dot. The dot charge couples to the elec- 
tric field associated with the bias voltage and leads to an 
electrostatic force (see Fig. [TJi) that acts on the movable 
dot with an electromechanical "shuttle instability"— as 
a possible consequence. In a magnetic shuttle device the 
dot magnetization (spin) in addition couples to the mag- 
netization of any nearby magnetic lead (or gate). This 
coupling gives rise to a short range exchange force (see 
Fig. Hb), which causes a spintronically induced mechan- 
ical response of the dot and hence one may talk about 
the spintro-mechanics of a magnetic shuttle device. 

In this Letter we will focus on one particular spin- 
tromechanical effect, viz. the formation of what we shall 
call "spin-polaronic states" in a magnetic shuttle device. 
In the limit of strong spintro-mechanical coupling the re- 
sult of a tunneling event that changes the net spin on 
the dot is a polaronic modification of the mechanical vi- 
brational states of the dot. The effect is similar to the 
modification caused by the addition of an electric charge 
on the do1ji£~— . In the latter case there is a Franck- 
Condon shift of the dot's position-^, while in the for- 
mer case the dot is displaced by the exchange force that 
appears as a spin is added in a combination of events 
that can be viewed as the formation of a spin-polaronic 
state. Since the tunneling matrix element is exponen- 
tially sensitive to the position of the dot one expects a 
spin-dependent exponential rcnormalization of the tun- 
neling probability caused by the formation of these spin- 
polaronic states. This is the origin of the exponentially 
strong spin-dependent tunneling effect which we will be 
discussed in detail below. 

The described spintro-mechanical effect on the shuttle 
device can be tuned by changing the bias voltage in order 
to inject n extra electrons onto the dot, thereby changing 
its spin as well as its charge. The diagram in Fig.[2]shows 
how the resulting shift of the equilibrium position of the 
movable dot (due to the exchange force) depends on its 
spin state. If the bias voltage is increased in order to go 
from the n = to the n — 1 state by adding one electron 
to the dot, the dot position will shift in different direc- 
tions with respect to the leads depending on whether the 
spin of this electron is up or down. This spin-dependent 
shift lifts the spin-degeneracy of the probability for tun- 
neling, which will be exponentially different for electrons 
with opposite spins. To see the effect of this spin discrim- 
ination of tunneling on the current, consider a magnetic 
shuttle device where the movable dot is closer to one of 
the leads (as indicated in Fig.QJi), which therefore is more 
strongly exchange-coupled to the dot than the other lead. 
If an electron with spin anti-parallel to the lead magneti- 
zation tunnels onto the dot the exchange-force will repel 
the dot from the nearby lead and decrease its distance 
from the more distant lead, which in turn will enhance 
the probability for this electron to tunnel to that lead and 
contribute to the current. If an electron with parallel spin 
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FIG. 2: Diagram showing how the equilibrium position of the 
movable dot in Fig. [Tb depends on its net charge and spin. 
The difference in spatial displacement discriminates transport 
through a singie-occupied dot with respect to spin. 

tunnels, the effect is the opposite and the probability for 
onward tunneling is decreased. The resulting current will 
therefore be spin polarized. 

If the bias voltage is further increased to add another 
electron to the dot, the resulting n = 2 state must be a 
spin singlet since we assume a single-level dot. This spin- 
zero n = 2 state would obviously not experience any ex- 
change force. However, in the Coulomb blockade regime 
the n = 2 state is energetically unfavourable, which is 
why (for an appropriate bias voltage) it can be expected 
to be frozen out even at relatively high temperatures. 
Spin-selective tunneling will therefore survive under the 
condition that Coulomb blockade of spintro-mechanical 
tunneling is preserved. 

So far we have argued that the injection of an extra 
electron onto the dot changes the exchange forced on 
the dot with the result that its equilibrium position is 
shifted by a distance A (see Fig. [5]). We will refer to the 
extra electron plus the shift of the dot position as a spin 
polaron. The Hamiltonian that describes the magnetic 
nanomechanical SET device has the standard form, ex- 
cept for its spin-dependent part (representing the mag- 
netic exchange energy) which now depends on the me- 
chanical displacement of the dot. Hence 

H = -ffloads + -^tunnel + -Hdot , (1) 

where 

-fflcads = ^ ajLcrakscrekscr (2) 

k,cr,s 

describes electrons (labeled by wave vector k and spin 
a =tj4-) m the two leads (s = L, R). Electron tunneling 
between the leads and the dot is modeled as 

■Htunnel = ^ T s (x) + B.C., (3) 

k,<r,s 

where the matrix elements T s (x) — exp(=Fx/A), with 
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A the characteristic tunneling length, depend on the dot 
position x. 

The movable single-level dot is modeled as a harmonic 
oscillator of angular frequency ojo, 

H dot = hujotfb + J2 n <? [ £ o - sign(er) J (x)] + Ecn^n^ , 

a 

(4) 

where sign(f, I) = ±1, Ec is the Coulomb energy associ- 
ated with double occupancy of the dot and the eigenval- 
ues of the electron number operators n a is or 1. The 
position dependent magnitude J(x) of the spin depen- 
dent shift of the electronic energy level on the dot is due 
to the exchange interaction with the magnetic leads (and 
any external magnetic field). 

The strength of the polaronic coupling a = A/x is de- 
termined by the ratio between the polaronic shift and the 
amplitude xq of the dot's quantum-mechanical zero-point 
oscillations. Another important parameter, /3 = A/A, is 
a measure of the effect of the spin-polaron formation on 
the probability for electron tunneling. As we will see later 
large values of both these parameters lead to a highly spin 
polarized tunneling current. 

If the polaronic shift A is small compared to the char- 
acteristic length I over which the exchange interaction 
changes significantly we may expand J{x) to linear order 
in x so that 



J (x) = J (0) + j x 



(5) 



and without loss of generality furthermore assume that 

„2 , 



= 0. In this case A = \j\xQ/(%u)o) and hence 



and P = a 



(6) 



Hluq A 

A full solution of the problem at hand can be obtained 
by solving the Liouvillc equation for the density matrix 
for both the electronic and vibronic subsystems. Two 
different limits determine different scenarios for the na- 
noelectromechanical response of the device. In the limit 
of low mechanical dissipation (high Q-factor) energy sup- 
plied by the external battery (used to maintain the bias 
voltage) may accumulate in the mechanical subsystem 
and eventually lead to an electromechanical instability 
and the onset of shuttle vibrations^. In the opposite 
limit of strong mechanical dissipation (low Q-factor), 
which we will consider here, the vibronic subsystem is 
kept in equilibrium at the ambient temperature T. 

In the weak tunneling limit a kinetic description of the 
electronic subsystem in terms of the probabilities P 7 s to 
occupy the on-dot electronic states is possible (7 and 
S are the eigenvalues of n-j- and respectively). By 
assuming that P 7 g only changes through single-electron 
tunneling events between the dot and one or the other of 
the two leads, one arrives at the rate equation 
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Here the tunneling rates T are are proportional to the 
squared modulus of the tunnel matrix elements between 
different stationary vibrational and electronic states. We 
are interested in the DC current through the device, 
which is determined by the time-independent solution of 
the kinetic equation ([7|). The spin-up current, e.g., can 
be expressed as 
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(8) 

Although a general solution of Eqs. (|7|) and (J5J for the 
stationary spin-up current is possible, the resulting ex- 
pression is quite complicated and not very suitable for 
a qualitative analysis. This is why we will here restrict 
ourselves to the experimentally most likely case of an 
asymmetric tunneling device with T^ s <C T?g. , s , . In 
this limit, where the electronic systems on the dot and 
in the nearest (right) lead will be in thermal equilibrium, 
the solution for the occupation probabilities P-ys will be 
independent of both T L and T R , while the current will 
depend on the tunneling rate T L to the most remote (left) 
lead. A direct calculation (following Ref. [l3h of the tun- 
nel matrix element rf . g, for example, gives the result 
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Here e = 6q — Ep/2, with Ep = a 2 ?iwo, is the po- 

laronically modified energy level on the dot, Wl = 
1 



T 



(0) 



with vl the electron states density at Fermi 

level in the left lead, is the tunneling-induced width of 
the dot level. Furthermore, 
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where /„(...) is a Bessel function, N is the average num- 
ber of excited dot- vibration quanta, 



N = 1/ [cx P (^o; /A;bT) - 1] , 



(10) 



and 



# = 1/ [cxp([e + sign(a)n^ - ^}/k B T) + 1] 



7<5;7'5 



Equation (|5J| for the current with the tunneling rates 
given by expressions like (|9]) and the occupation 
probabilities Ps-y obtained from the kinetic equation ([7]) 
gives an analytical expression for the current. At low 
temperatures and low bias voltages the current is affected 
by both the familiar (spin-independent) polaronic block- 
ade of tunneling (through the parameter a)^2 ~ 15 ' 17 and 
by spin-selective tunneling (through /3). For example, at 
V ~ and T « the partial current of electrons with 
spin a is 



e — exp U 



sign(a)/3 



(11) 
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where a spin-independent prefactor of order 1 has been 
omitted. 

An increase of voltage or temperature lifts the po- 
laronic blockade of tunneling when max{eVy2, fc^T} > 
Ep. Tunneling will still be spin selective, however, as 
long as max{eV/2, kpT} < Ep + Ec so that double oc- 
cupancy of the dot is prevented by the Coulomb blockade. 
In this case the expression for the spin currents take the 
form 

h ~ exp ([2N + 1] |f - 2 sign(a)^) , (12) 

where N is given by Eq. (fTU|) . At even higher volt- 
ages or temperatures, max{eV/2 7 fcsT} > Ep + Ec, the 
Coulomb blockade of double dot occupancy is also lifted. 
A doubly occupied dot level will have no net spin and 
therefore the exchange force on the dot vanishes, no spin- 
polaron is formed, the spin currents are independent of 
spin projection, ft — in Eq. (|12[) . and the spin-polaronic 
stimulation of the current ends with the result that the 
current drops at eV ~ 2(Ep+Ec) (at low temperatures). 
For higher voltages the current will be oc exp(2T/T*), 
where k B T* = Hujq/xq (assuming k B T > Tiujq). This is 
because the amplitude of the random thermal vibrations, 
which effectively reduce the tunneling distance, increases 
with temperature. 

The higher current at intermediate bias voltages is an 
example of "Coulomb promotion of tunneling" , which 
here depends on the spin-dependent shift of the dot po- 
sition associated with the formation of spin polarons. 
Coulomb promotion of tunneling may also appear due 
to the standard spin dependent tunneling phenomenon 
(mentioned in the introduction) as discussed in Ref. fill 

For a numerical example we consider an extended 
quantum dot in the form of a single- wall carbon nanotubc 
suspended between two electrodes (source) and probed 
by a magnetic-metal STM tip (drain), which is similar 
to the set-up used by LeRoy et al*&. For a suspended 
tube length of order 1 yum the vibration frequency ujq /2tt 
of the fundamental bending mode is of order 100 MHz 
and hence Hluq ~ 1 /zeV (Tiwo/fcs ~ 0.01 K). Further- 
more, assuming a diameter of 2 nm, the suspended mass 
M of such a tube is such that the quantum oscillation 
amplitude xq ~ y / h,/2Mu>o ~ 0.01 nm. The tunneling 
length is an atomic distance, so A ~ 0.1 nm is a reason- 
able estimate and hence A ~ 10 .To S> xq. Finally, if one 
approximates the spatial gradient of the exchange field 
j = d.J(x)/dx ~ J/A where J - 0.1 meV (J/k B ~ 1 K) 
one finds that a ~ 10 and ft ~ 1. These parameter values 
are large enough for the predicted effect to be strong as 
is evident from Fig. |3l where the current- voltage charac- 
teristics and spin polarization of the current are plotted 



as functions of voltage and temperature. In particular 
it is clear from Fig. [3] that the spin polarization of the 
current can be very close to 100%. The characteristic 
temperature T* is a few kelvin with these parameters. 

One notes that the charging energy for the extended 
nanotubc quantum dot of our example is rather small, 
4 K, and hence well below room temperature. However, 
one may consider using functionalized nanotubes^ or 
graphene ribbons^ with one or more nm-sized metal or 
semiconductor nanocrystals attached. The small size of 
the crystals could give several orders of magnitude larger 
charging energies while not much affecting the low me- 
chanical vibration frequencies of the carbon resonators. 
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FIG. 3: Spin polarization of the current through our model 
magnetic NEM-SET device for hu> /k B = 0.01 K, A = 
10 xo ~ 0.1 nm, Aj'/fc s = 2 K [hence E P = 400 fiuo] and 
Ec = iOOhiuo: (a) ln(Ii/If) as a function of bias voltage 
V and temperature T (b) ln(/.|.//-|-) as a function of V at 
T — 30 hujo/kB [grey (red) curve] and 1 = 1 (V) for the same 
parameters, i.e. along the dashed line in panel (a) [black 
curve] . 

In conclusion wc have demonstrated that a 
spintro-mechanical coupling caused by spatially 
non-homogeneous exchange interactions in mag- 
netic nanoelectromechanical single-electron tunneling 
structures may result in a very strong spin dependence 
of the probability for electron tunneling. As a result 
the theoretical value of the spin polarization of the 
electrical current can essentially be 100%. For future 
work we note that the strong spin-dependence of me- 
chanical displacements demonstrated here in the limit 
of strong mechanical dissipation, may lead to nontrivial 
spin-mcchanical dynamics in the weak dissipation limit 
if both the spin dynamics and the mechanics is coherent. 
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